Is an Ultra-Cold Strongly Interacting Fermi Gas a Perfect Fluid? 



J. E. Thomas 

Physics Department, Duke University, 
Durham, NC 27708-0305, USA 



Abstract 

Fermi gases with magnetically tunable interactions provide a clean and controllable laboratory 
system for modeling interparticle interactions between fermions in nature. The s-wave scattering 
length, which is dominant a low temperature, is made to diverge by tuning near a collisional 
(Feshbach) resonance. In this regime, two-component Fermi gases are stable and strongly in- 
teracting, enabling tests of nonperturbative many-body theories in a variety of disciplines, from 
high temperature superconductors to neutron matter and quark-gluon plasmas. We have devel- 
oped model-independent methods for measuring the entropy and energy of this model system, 
providing a benchmark for calculations of the thermodynamics. Our experiments on the expan- 
sion of rotating strongly interacting Fermi gases in the normal fluid regime reveal extremely low 
viscosity hydrodynamics. Combining the thermodynamic and hydrodynamic measurements en- 
ables an estimate of the ratio of the shear viscosity to the entropy density. A strongly interacting 
Fermi gas in the normal fluid regime is found to be a nearly perfect fluid, where the ratio of the 
viscosity to the entropy density is close to a universal minimum that has been conjectured by 
string theory methods. 



1. Introduction 

Tabletop experiments with degenerate atomic Fermi gases near a Feshbach resonance |ll 2] 
provide models for strongly interacting Fermi systems in nature. Feshbach resonances yl, \4\ 
arise when different hyperfine channels have different magnetic moments, as occurs in 6 Li and 
40 K atomic Fermi gases. A bias magnetic field is applied to tune the total collision energy of 
the incoming continuum state into resonance with that of a bound state in a closed channel. At 
resonance, the zero energy s-wave scattering length as diverges, and the collision cross section 
is only limited by unitarity, i.e., cr oc A 2 B , where Xb is the de Broglie wavelength. Even though 
it is a dilute system, a unitary atomic gas is the most strongly interacting non-relativistic system 
known |5fl. 

Strongly interacting Fermi gases exhibit strong pairing interactions, of interest in the field 
of high temperature superconductivity JfJ, neutron stars, and nuclear matter 10, B S OH- The 
common feature which all of these systems share is a strong interaction between pairs of spin- 
up and spin-down fermions. The strongly collisional normal fluid exhibits extremely low vis- 
cosity hydrodynamics and elliptic flow 11Q, analogous to the hydrodynamics of a quark-gluon 
plasma lUl ll2H . In contrast to other Fermi systems, atomic gases enable magnetically tunable 
interactions [13, 14, lit], variable energy ['3 12, 18], and variable spin populations 



Fermi gases near Feshbach resonances are now being widely studied ['21]. 
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It has been suggested that a strongly interacting Fermi gas may serve as a model for the most 
"perfect" fluid, which has a minimum viscosity. A simple argument based on quantum mechanics 
places a lower bound on the viscosity [ 22] . The shear viscosity is of order n oc n pA m f p , where n 
is the density, p is the average momentum and A m f p is the mean free path Il23n . The Heisenberg 
uncertainty principle requires pA m f p > H, so that the shear viscosity satisfies n > nh. In high 
energy physics, where particle number is not conserved, the ratio of the viscosity to the entropy 
density s ^ nk B is often considered. Hence, one expects that n/s > h/k B - 

Recently, using string theory methods, it has been conjectured that there is a universal strong 
coupling lower bound W2411 . 

n ->^. a) 

s An k B 



Although one can imagine high entropy systems which might violate the lower bound 12511 . cur- 
rently no fluid that even achieves the lower bound is known. If the viscosity conjecture is correct, 
it represents an extremely important advance in the understanding of many-body physics 112511 . 
Hence, it is of great interest to explore minimum viscosity quantum hydrodynamics in strongly 
interacting Fermi gases as a model system [5,2l. 

An important feature of strongly interacting Fermi gases is the property of universality H, 



211 128L 12911 . The system exhibits scale invariance, in the sense that, at zero temperature, the 
interparticle spacing L sets the only microscopic length scale at resonance, leading to universal 
behavior. In a uniform strongly interacting gas, the ground-state energy is a universal fraction, 
denoted 1 + J3, of the energy of a noninteracting gas at the same density Illl30|]. This universal 
energy relationship was originally explored theoretically in the context of nuclear matter [0, [gl 
BO and has now been measured using ultracold Fermi atoms El BE El [II. Our best 
current measurements of J3 are given in Ref. 113 311 . 

Universality automatically leads to the natural "quantum viscosity scale" hn ll34ll . The shear 
viscosity has natural units of momentum/area. In a strongly interacting Fermi gas, the natural 
momentum is h/L, while, the natural area is the unitary collision cross section, of order L 2 for 
temperatures at or below the Fermi temperature. The shear viscosity is then of order h/L 3 or hn. 
It is therefore natural to write the viscosity in the form, 

n = ahn, (2) 

where a is a dimensionless parameter, which is generally a function of the local reduced tem- 
perature, T/T F (n), where T F (n) is the local Fermi temperature, oc n 213 . We note that for water 
a 300, while for air, a 6000. For liquid He near the /1-point, a 1, in the quantum regime. 
As we will see, our estimates for a strongly interacting Fermi gas are significantly lower. 
For a trapped gas, we are able to estimate the ratio n/s as 112611 

n = h (a) ^ 
s k B Slk B ' 

where S /ks is the average entropy per particle for the trapped gas and (a) is the trap averaged 
shear viscosity in units of hn. In the following, we describe our recent studies of the thermody- 
namics, i.e., the measurement of the energy and entropy. Then we will describe our studies of 
the hydrodynamics and the estimate of the shear viscosity. Using these results, we compare to 
the lower bound of Eq. [T] 
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2. Measuring the Entropy and Energy 



Recently, we have developed model independent methods for measuring the entropy and 
energy of a strongly interacting Fermi gas il8ll22l]- Energy measurement is based on the virial 



theorem 12811 . Since the interparticle spacing and the thermal de Broglie wavelength are the 
only length scales when the cloud is tuned to a broad Feshbach resonance, the local pressure 
is a function only of the local density and temperature. In this case, one easily verifies that 
the virial theorem holds using elementary thermodynamic arguments, which is confirmed by 
experiment 1 28] . The atoms are confined in an optical trap, which at low temperatures, provides 
a nearly harmonic trapping potential, yielding the energy 

E = 2(U) = 3moj 2 z (z 2 ) s . (4) 

Here, we have assumed a scalar pressure, so that the harmonic trapping potential energy is iden- 
tical in all three directions, x, y, z. We make measurements of the mean square size in the long 
z-direction of the cigar- shaped cloud, Fig. [T] The spring constant, moj 2 is typically determined 
within 0.5%, and the mean square size of the strongly interacting gas (z 2 )s is determined within 
2%. Despite the fact that the gas generally contains condensed superfluid pairs, non-condensed 



Figure 1: Absorption image of an optically trapped, strongly interacting Fermi gas. The length of the cloud is ^ 200 /urn. 

pairs, and unpaired atoms, all strongly interacting, this remarkable result shows that a simple 
measurement of the mean square size of the trapped cloud determines the total energy. 

The entropy is measured by means of an adiabatic sweep of the bias magnetic field, to tune 
the scattering length as from the strongly interacting regime, na\ » 1, to the weakly interacting 
regime, where na\ « 1. In the weakly interacting regime, the entropy Sw of the cloud is nearly 
that of an ideal Fermi gas in a harmonic trap, where 

S w ^S ideal = J deD(e)s(e,T). (5) 

Here s(e, T) is the Boltzmann entropy for an orbital of single particle energy e at temperature T, 
and D(e) is the density of states for the trap. The chemical potential is determined by normalizing 
the total occupation number to the total number of atoms. Using this, the entropy, energy and 
mean square size of the cloud are then determined for a given temperature T. Eliminating the 
temperature, the entropy is given in terms of the mean square cloud size of the weakly interacting 
gas (z 2 )w, which is readily measured, 

S\y - S ideal ((Z 2 )w ~ (z 2 )wo) • (6) 
3 



Here, (z 2 )wo is the mean square size of the cloud in the ground state, which is estimated by 
fitting a Thomas-Fermi density profile to the density profiles at the lowest temperature, yielding 
the Fermi radius, cr z , from which (z 2 )wo = The measured value is in very good agreement 

with calculations based on the local chemical potential at zero temperature fio(n) ll33ll . We note 
that writing the entropy as a function of the mean square size relative to the ground state assures 
Sw = for the measured ground state cloud size. This method also improves the ideal gas 
approximation by suppressing small mean field corrections to the ideal gas cloud sizes at finite 
interaction strength. We find that the corrections to the ideal gas entropy for the finite interaction 
strength are within a few percent except at our lowest temperatures, where the correction is 
^ 10% Q. 

To verify that the sweep of the bias magnetic field is adiabatic, we note that after a round trip 
sweep lasting 2 s between the strongly and weakly interacting regimes, the energy is found to be 
within 2% of that obtained by simply holding the strongly interacting cloud for 2 s. Hence, 

Ss=S w . (7) 

To perform the measurements, an atom cloud is first cooled by lowering the trap depth to 
achieve forced evaporation to a temperature near the ground state. Then energy is added by 
releasing the cloud for a precisely controlled time and then recapturing the cloud. This method 
reproducibly adds energy to the cloud, which is allowed to equilibrate. The energy of the strongly 
interacting cloud is then measured from the cloud size. A second cloud is then created using the 
same parameters as the first. The bias magnetic field is swept over 1 s to the weakly interacting 
regime, where the mean square size is measured. Together, these measurements yield the energy 
and entropy of the strongly interacting gas, Fig. [21 

The scaling of the energy with entropy is quite different for low energies E/Ep < 0.8 than 
for higher energies E/E F > 0.8. We attribute this change in the thermodynamics to a superfluid 
transition. Fitting a smooth curve to the data, and using T = dE/dS , we find that superfluid- 
normal fluid transition occurs at a temperature T - 0.2 7> 113 311 . which is quite high. In a charged 
condensed matter system, where 7> corresponds to an eV or 10 4 K, this would correspond to a 
superconductor transition at 2000 K! 



3. Estimating the Shear Viscosity 

For a unitary Fermi gas, the bulk viscosity is believed to vanish J35, 36], so that the shear 
viscosity determines the hydrodynamic damping rate. Our first estimates of the shear viscosity 
were made by measuring the damping rate of the radial breathing mode 12611 . Schafer ||5|] used 
this damping data to estimate the shear viscosity and combined the results with our entropy data 
to make a comparison with the lower bound of Eq.Q] We noted previously that edge effects in the 
trap, such as interactions between the hydrodynamic and ballistic components of the cloud, might 
have increased the observed damping rate in the trapped gas. The damping rates do not appear to 
scale properly with atom number at fixed E/E F if viscosity is assumed to be the primary cause 
of damping II26I1 . However, the estimated viscosity is in the quantum regime. 

Recently, we have estimated the shear viscosity in a different way, by measuring the expan- 
sion dynamics of a rotating Fermi gas, which is released from the optical trap 113 711 . In this case, 
the gas expands freely, eliminating the direct effects of the edges of a trapped cloud. The gas 
is cooled by evaporation to near the ground state and a controlled amount of energy is added. 
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Figure 2: Measured total energy per particle in units of Ep of a strongly interacting Fermi gas at 840 G versus its entropy 
per particle in units of For comparison, the dot-dashed green curve shows E(S) for an ideal Fermi gas. Ep is the 
Fermi energy of an ideal Fermi gas at the trap center. 



Then the trap is rotated abruptly to excite a scissors mode. The gas is released and imaged af- 
ter a selected expansion time. Fig. [3] shows typical data, for different initial angular velocities. 
The angle of the long principal axis of the cloud is measured as a function of time after release. 
The data reveal that as the gas expands, the angular velocity increases, which is a consequence 
of irrotational hydrodynamics: The moment of inertia decreases as the aspect ratio approaches 
unity. 

Remarkably, the cloud for the normal fluid at E/E F = 2.1 behaves almost identically to the 
superfluid cloud for E/Ep = 0.56. Indeed, the effective moment of inertia is quenched well 
below the rigid body value in both cases, and is in very good agreement with expectations for 
irrotational flow ll37ll . 

Irrotational flow is expected for the superfluid, since the velocity field is the gradient of the 
phase of a macroscopic wavefunction. However, irrotational flow in the normal fluid requires 
very low shear viscosity. As energy is added to the gas, we find that the expansion dynamics 
slows down compared to ideal, isoentropic irrotational flow. 

To estimate the shear viscosity, we use a simple model. We assume that the slowing down 
of the dynamics, compared to ideal irrotational flow, arises from shear viscosity. In the absence 
of viscosity, the expansion is isoentropic. In this case, for release from a harmonic trap, an exact 
solution to the hydrodynamic equations is obtained for a velocity field that is linear in the spatial 
coordinates ff3vh . We then add to the hydrodynamic equations a term that is the divergence 
of the pressure tensor arising from shear viscosity [ 38] . The time evolution depends on the 
trap average (a), Eq. [2 which increases as the energy is increased. At energies just above the 
superfluid transition, the viscosity is sufficiently small that the uncertainty is determined by the 
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Figure 3: Hydrodynamic expansion of a rotating, strongly interacting Fermi gas, for different initial angular velocities. 

accuracy of the trap parameters (oscillation frequencies in the three directions are determined 
within 0.5%). Note that for the lowest viscosities, the fit can return a negative value of (a), 
which is an artifact arising from measured trap parameters that predict a slightly slower evolution 
for perfect irrotational flow than that measured. Fig. [4] shows how the estimated shear viscosity 
depends on the energy of the cloud. The shear viscosity is given in units of the quantum viscosity, 
i.e., in units of hn, where n is the density. The red (blue) data are taken at trap depths that are 
20% (5%) of the maximum attainable. The agreement in the data for both trap depths shows that 
anharmonicity in the Gaussian trapping potential is not significant. 
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Figure 4: Estimated shear viscosity in units of hn versus energy in units of Ep. Red (blue) circles denote data at 20% 
(5%) of maximum trap depth. 



By combining the entropy measurements with the viscosity estimates, we are able to estimate 
the ratio of the shear viscosity to entropy density. Fig. \5\ shows how the results compare to the 



string theory conjecture for the minimum ratio, Eq. [TJ Our estimates of the viscosity suggest 
that a strongly interacting Fermi gas in the normal fluid regime (above 0.8 E F ) is a nearly perfect 
fluid. 
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Figure 5: Estimated ratio of the shear viscosity to the entropy density. Dotted line shows the string theory conjecture [24] 
for the minimum ratio, Eq.^ 
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